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ABSTRACT 

The Levy stability analysis is carried out for e+e~ collisions at Z° 
mass using Monte Carlo method. The Levy index n is found to be 
H = 1.701 ± 0.043. The self-similar generalized dimensions D{q) and 
multi-fractal spectrum /(a) are presented. The Renyi dimension D{q) 
decreases with increasing q. The self-similar multi-fractal spectrum is 
a convex curve with a maximum at g = 0, a = 1.169 ± 0.011. The 
right-side part of the spectrum, corresponding to negative values of q 
is obtained through analytical continuation. 
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1 Introduction 



The study of fluctuations in particle physics already has a long history going back to early 
cosmic-ray observations [Q]. Bialas and Peschanski proposed to study non-statistical fluctu- 
ations in multi-particle production by the method of factorial moments 

where a region A in 1-, 2- or 3-dimensional phase space is divided into M cells, rim is the 
multiplicity in the mth cell, and (■ ■ ■) denotes vertically averaging over the event sample. 
There are a large variety of experiments showing that Fq{M) possess an anomalous scaling 
behavior |Q 

Fg{M) oc (M)*(«) (M-^oo) , (2) 

Note that in e^e^ collisions the fluctuations or power-law scaling exist in 3-dimensional 
phase space and the dynamical fluctuations has been proved to be isotropic [QPP. This 
means that the multi-hadron final states in these collisions are self-similar fractal. 

Since the fractal property of the multi-hadron final states in e"'"e~ collisions has been 
established, the next step is to study this fractal system in more detail. The aim of the 
present paper is to perform the Levy stability and multi-fractal spectrum analysis for the 
fractal system of the multi-hadron final states in e"'"e~ — collisions at the Z° energy — 
y/s = 91.2 GeV using Monte Carlo method. 

2 Measurement of intemittency (l){q) 

In our Monte Carlo simulation a total number of 200000 events of e"'"e~ collisions at 91.2 
GeV are produced by JETSET 7.4 generator. The primary produced partons considered is 
a quark-antiquark pair, moving back to back in the rest frame. This implies a cylindrical 
symmetry about the quark direction. An appropriate frame to study the development of the 
qq system is therefore the thrust frame 0, the appropriate variables are {y,pt,ip). 

In order to avoid the influence of a non-flat distribution of the variables y, pt, and (p 
on the investigation of the dynamical fluctuations, all variables are transformed into their 
corresponding cumulant forms [§. 

ly' P{y)dy jpll P{Pt)dpt Q PW)d'p 

The results for the three dimensional FM defined in Eq.(l) are shown in Fig.l. After 
omitting the first point to eliminate the influence of momentum conservation [10|, the 



three dimensional 2 ~ 7 order InFg's versus InM data are fitted very well by a straight line 

InF, = 0(g)lnM + 6, (4) 



2 



where (/)(g) is the intemittency index. The resulting parameters and the corresponding x^/DF 
are hsted in table I. 

Table I The fit parameters of 3-D FM for the Jetset7.4 



order q 


0(9) 


b 




2 


0.191 ±0.001 


-0.304 ±0.003 


57/9 


3 


0.588 ±0.003 


-0.914 ±0.009 


27/9 


4 


1.135 ±0.009 


-1.690 ±0.023 


20/9 


5 


1.784 ±0.024 


-2.554 ±0.061 


7/9 


6 


2.435 ±0.091 


-3.289 ±0.237 


4/8 


7 


3.223 ±0.197 


-4.221 ±0.519 


3/4 



The results show power-law scaling when phase space is partitioned isotropically, indicat- 
ing the existence of self-similar (isotropic) dynamical fiuctuations. Thus the multi-hadron 
final states in e"'"e~ collisions at y/s = 91.2 GeV is a self-similar fractal. 
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Fig.l The three-dimensional 2 ~ 7 order factorial moments as functions of the partition number M 



3 Levy index 

As a self-similar fractal system, the multi-particle final states in high energy collisions can be 
charactered by an important parameter — the Levy stability index /i. This parameter tells 
us the behavior of elementary fiuctuations at the tail. It is also known as the degree of multi- 
fractality, /i = for mono-fractals, /i < 1 corresponds to the so-cahed "calm" singularities, 
while /i > 1 correspond to the "wild" singularities. 

The intermittency index 0(g) is related to the anomalous dimension dq as 

<P{q) = {q- l)dq, (5) 
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and the later is expressible by the Levy index /i as 



da 

" fi-lq-1 
or in terms of the intermittency index 

0(g) = 0(2) 



(Ci = const), 



- 2' 



(6) 



(7) 
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Fig. 2 The intermittency index 0g as a function of order q. 
SoUd curve is from Eq.(9), dashed curve is from Eq.(8). 

According to the Levy stability theory, the Levy index /i is confined to < ;U < 2, and 
the Central Limiting limit corresponds to ;U 2. In this limit Eq.(7) passes to 



0(g) = 0(2) 



g(g- 1) 



(8) 



Note that Eq.(8) is valid for any real value of g, including g < 0, while Eq.(7) in the present 
form is inapplicable to g < 0. An analytical continuation of Eq.(7), which is valid for g < 
and will, in the limit /i 2, pass to Eq.(8) in the whole region of g is 



0(g) = 0(2) 



(g^)^/^ - g 
2^* - 2 



(9) 



The distribution of 0(g) versus g is shown in Table. I and Fig 2. We fit it to Eq.(7), 
shown as solid curve in the figure. This curve has been extrapolated to g < 1 using Eq.(9). 
Also shown in the figure is the result from Eq.(8) (dashed curve), which corresponds to 
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the Central Limiting case. Note that the g < 1 part is needed for getting the multi-fractal 
spectrum to the right of its maximum, cf. Fig.3. 

The Levy index fi for the hadronic system in e~^e~ collision at y/s = 91.2 GeV obtained 
from the fit is 

/i = 1.701 ± 0.043. (10) 

This value of /i lies in the region [0, 2] in consistent with Levy stability. Since /z is greater 
than unity, the multi- fractal corresponds to "wild" singularities. 



4 Multi- fractal spectrum and Renyi dimensions 

The multi-fractal method is a widely used tool in many branches of physics and science in 
general |1T[. The Renyi dimension Dq {= 1 — dq) and multi-fractal spectrum f{a) are often 
used to study the multi-fractal. The D{q) (also denoted as Dq) depends on order q, being 
a decreasing function of q for multi-fractals in general. The Renyi dimension is sometimes 
called as generalized dimension. The multi-fractal spectral function f{a) and the Renyi 
dimension (generalized dimension) D{q) for hadronic system are related to intermittency 



index 0(g) as [12 



r(g) = g - 1 - 0(g), a(g) = 

dq 

/(a) = ga-r(g), Diq) = (11) 

Where D(q)\q=Q = Dq = Dp is the fractal dimension, D{q)\q=i = Dj is the information 
dimension and D{q)\q=2 = D2 = a;(2) — f{a{2)) = z/ is the correlation dimension. 

The continuous order index 0(g) are calculated by Eq.(9). Then the Renyi dimension 
D{q) and multi-fractal spectrum /(a) are calculated using Eq.(ll) with g = -1 to 6.8, step 
0.2. The -D(g) versus g and the f{a) versus a are shown in fig.3(a) and (b). 
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Fig. 3 (a) The Renyi dimension D{q) as a function of order q 
(b) The multi-fractal spectrum f(a) as a function of a 

It can be seen from Fig.3 that, the Renyi dimension D{q) decreases from 1.17 to 0.437 
as the order q increases from —1 to 6.8 which means that there is a multi-fractal behavior in 
multi-production at e+e~ coUisions. The multi-fractal spectrum f{a) is lying between —2.10 
and 1 in the above-mentioned a region, having a maximum at q — 0, a — 1.169±0.011. The 
corresponding fractal dimension is Dp = f{a{Q)) = D{0) = 1.003 ± 0.001. 

The tangents to the curve /(a) = a with slopes d/(Q;)/da; = 1 and 2 determine the 
points q{1) = 0.882 ± 0.006 and «(2) = 0.701 ± 0.006, indicating that the information 
dimension of the system is Dj = D{1) = 0.882 ± 0.006, and the correlation dimension is p 
= 2a(2) - /(a(2)) = D{2) = 0.787 ± 0.012. 

The Renyi dimensions D{q) could also be directly calculated using Eq.(ll). The resulting 
values are listed in Table II together with those obtained from the Multi-fractal spectrum 
as described above. It can be seen that the consistency is good. 

Table II The Renyi dimensions obtained from two different methods 



method 


Dp{D{0)) 


Di{D{l)) 


u,{D{2)) 


Equation (11) 


1.000 ±0.000 


0.881 ± 0.027 


0.788 ± 0.003 


Multi-fracltal spoctruni 


1.003 ±0.001 


0.882 ± 0.006 


0.787 ± 0.012 



6 



5 Conclusions 



A self-similar multi-fractal analysis using MC method in e+e^ collisions at ^/s = 91.2 GeV 
is performed basing on the measurement of factorial moments. The Levy stability analysis 
is carried out and the Levy index in found to be = 1.701 ± 0.043 consistent with Levy 
stability. The multi-fractal spectrum f{a) are presented which is a convex curve with a 
maximum ai q = 0., a = 1.169 ± 0.011. The Renyi dimensions D{q) are obtained both from 
their relation with the intermittency indies 4>{q) and from the multi-fractal spectrum. The 
results are in good consistence. 
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